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Abstract 

The two-point Green function of the massive scalar (3 + l)-quantum held the¬ 
ory with A ip 4 interaction at finite temperature is evaluated up to the 2nd order of 
perturbation theory. The everaging on the vacuum fluctuations is separated from 
the everaging on the thermal fluctuations explicitly. As a result, the temperature 
dependent part of the propagator is expressed through the scattering amplitudes. 
The obtained expression is generalized for higher orders of perturbation theory. 
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1 Introduction 


In 1969 R. Dashen, S. Ma and H.J. Bernstein have suggested [1] a generalization of the 
Beth-Uhlenbeck formula. They have obtained the complete virial expansion for the grand 
potential where the n-tli virial coefficient was expressed in terms of special traces of 
the n —> n ,9-matrix elements 

/ , (Si - Sio) = T£ »/ dEe -^( TlnAS -^s), (i) 

n ~^ n-m 

( T = 1/(3 is the temperature, z is the activity, S is the scattering matrix, A is the exchange 
operator, m is the mass of a particle and Ho is the grand potential of the free particle 
system). The compact and invariant form of the eq.(l) was the reason to claim [1] the 
validity of this representation also in the relativistic case. We share this point of view in 
spite of absence of any direct quantum field theory derivations. 

In 70-th we exploited this formalism - so called S-matrix formulation of statistical 
mechanics - for the phenomenological investigation of the hot hadron matter thermody¬ 
namics [2-5]. The equation of state derived on the basis of Regge phenomenology for 
the scattering amplitudes appeared to be essentially nonideal at GeVs temperatures. We 
have used it for some astrophysical and cosmological applications [6-8]. In particular 
due to a Van-der-Waals-like nonideality of our equation of state we have obtained in 
the framework of the Fridman model the exponential expansion of the Universe in the 
vicinity of the phase transition, providing for a solution the of horizon, flatness, isotropy, 
primary fluctuations and some others problems. A few years later Guth published [9] his 
outstanding Inflationary Universe scenario based on a quantum field theory approach. 

At present the general fashion for the theoretical treatment of hadron matter (quark- 
gluon plasma) is ingenuous QCD: perturbation theory, lattice simulations. Nevertheless 
we believe that the A-matrix approach to the problem is fruitful at least in a qualitative 
and heuristic sense (as for quantitative analysis see for example ref. [10]). To be more 
acceptable S'-matrix approach needs of course in detailed QFT tests. The perturbation 
theory analysis is a step in this direction. 

In this work we carry out two-loop calculations for the temperature dependent part of 
the 1PI two-point Green function of A tp 4 QFT - the simplest but nontrivial quantum field 
model. We find that it can be represented through the thermal everages of corresponding 
renormalized scattering amplitudes. 

2 The model 

We considered a scalar held (p(x) in a box with the volume V at the temperature T = 
(3~ 1 . The points of the Euclidian space-(imaginary)time are enumerated by coordinates 
x = (x 0 ,x). The scalar product is 

(xy) = x 0 y 0 + (xy), x 2 = x 2 0 + x 2 . 
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The field (p(x) is periodical along the time direction ip{xo,x) = <p(xo + /3,x) according 
to Bose-Einstein statistics. Space boundary conditions are irrelevant in thermodynamical 
limit V —> cx). 

The action is 


£[</?] = Jd 3 x J dxo(jCo(p) + £/(<d)) 


where the free-field Lagrangian is 


The interaction due to 


A)0) = ^(x)(m 2 0 - n)(p(x), 

d 2 d 2 d 2 d 2 

dxl dx 2 dx 2 dx 2 ' 




where mo, Aq are the bare mass and the coupling constant. The generating functional is 


exp{Z[J]} = J 


where J(x) is the external source and 


(J,<p) = j d 3 x j dx 0 J(xMx). 


The Green functions are defined as usual [11,12] 

G (») (ll * ) = __ 

The one particle irreducible (1PI) Green functions are 


Y [n \x x ,,.. : x n )- 


H(x 1 ). ..5ip(x n y 


where 


= sz M sjm = _ Jix) 

5J(x )’ 8'4>{x) 


We use the Fourier transformation (direct and inverse) in the appropriate normalization 

f(x„,x) = „,p), 

\ Z7r ) p n J 


f(Po,p) = / d 3 x 


p 0 = 2nlT, 1 = 0, ±1, ±2,.... 
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3 1PI r^-function 


At T = 0 and V —> oo the Green functions depend on momenta only. At T / 0 
they depend also on T. The T^-function is the inverse propagator: G^ ■ T^ 2) = 1. The 
observable corresponding to the propagator is the space (equal-time) two-point correlation 
function 


«(N) 


T „ r cPpe ipx 

O0 3 lW r W(p 2 0 ,p 2 ,T )' 


( 2 ) 


The correlation length x c (T) is determined by the nearest to the real axes zero of the 
denominator of integrand ( 2 ) in the complex plane | p \ at po = 0 . 


T (2) (0, -x;\T),T) = 0. 


At T = 0 the correlation length is simply the inverse mass of the boson and at T ^ 0 it 
defines the Debye screening. 

The goal of this paper is to separate the vacuum (quantum) fluctuations from thermal 
ones. Let us remind the well known one-loop result for T®-function 


F i2) {p 2 0 ,p 2 ,T) =p 2 + m 2 0 + 


A 


d 4 q 


+ 


A 


d 3 q 


2 (27t) 4 J q 2 + m 2 (2n) 3 J 2uj{ed u> — 1) ’ 


(3) 


where 

2 2 i ->2 2 —*2 i 2 

P = Po + P i U ~ P + m ■ 

This expression clarihes what we mean in “vacuum” and “thermal” fluctuations. We call 
the fist integral in (3) “vacuum loop” and the second one - “thermal loop”. The full loop 
is divided in two terms - additivcly in the first order of perturbation theory 




Thermal fluctuations exponentially vanish at T 0 and 

^ (2) (pI,p 2 ,t = 0 ) = T ( 2 ) ( p 2 ) 

where T^(p 2 ) is the usual inverse propagator. Let us define the temperature dependent 
part by extracting from the full propagator its value at zero temperature 

i(pI,p 2 ,t) = y^\pIp\t) - r ( 2 ) o » 2 ). ( 4 ) 

Since the 2 —> 2 scattering amplitude in the 1st order of PT is 

A (2 \p,q-, p',q') = -A 


the 7 -function can be written as 


7(Po,2,p 2 ,T) 


1 f d 3 qA < ' 2 ' ) (p,q; p,q) 

( 2 vr ) 3 J 2iu{el 3w - 1 ) 


(5) 


Of course, this expression seems artificial for the trivial 1st order of PT but it shows what 
we search for higher orders. 
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4 The second order of PT 


The calculations at T ^ 0 are more complicated in 2nd order compared to that at T — 0 
by two reasons: the regularization (cut off) is not symmetric in momentum space and 
there is a loop summation instead of loop integration. The last obstacle is important 
because at finite temperature there appears loop divergences depending on T. They come 
from different diagrams. One of them contains a simple sum and the other - the double 
sum. So care should be taken when transforming sums to integrals to cancel the above 
mentioned divergences. 

The T^-function in the 2nd order is saturated by following diagrams 



(7) 

The appropriate regularization is ] p \ < A, i.e. it is a cylinder in momentum space 
rather then a sphere as usual at T = 0. Loop summations over zero’s components of 
momenta are convergent by itself. Subtraction points are defined at T — 0 as follows 


r ( 2 V) 



d r< 2 V) 

dp 2 



( 8 ) 


T (4) (Pi,P2,P3,P4) 


= -A, 

Pi =o 


(9) 


where m and A are renormalized (observable), mass and coupling constant. 

The technical problem is to transform the loop sums into integrals. For the simple 
sum over go the answer is [12] 


Sv = T*£ 


Qo 


(g 2 + m 2 ) 1 


1 

27T 


dq o 


(g 2 + rri 2 ) 1 


1 + 


o 


,v=l,2, 


( 10 ) 


or 

J_ r dq 0 1 

1 2n J Q 2 + m 2 — 1) ’ 

S' = — [ d% 1 \ 

2 2tc J Q 2 + m 2 u dio ^^(e^ — l)y 

There is more trouble with double sum S :i contained in the sunrise diagram 


( 11 ) 

( 12 ) 


O 
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( 13 ) 


c = 7 12 V ___ 

3 ^ ( q 2 + m 2 )(k 2 + m 2 )[(g + k — p) 2 + m 2 } 

The substitution (10) is not correct in this case owing to the pole surface in the complex 
manifold go <S> fco caused by the factor [(g + k — p) 2 + m 2 ] -1 in (13). Careful manipulations 
give the result 

„ _ 1 f _ dq 0 dk 0 _ 

3 (2n) 2 J (g 2 + m 2 )(k 2 + m 2 )[(q + k — p) 2 + m 2 ] 

3 1 r d% ( 1 1 \ 

2oj k (eP Wk — 1 ) 2n J (q 2 + m 2 ) y (g + k — p) 2 + m 2 (q — k — p) 2 + m 2 J 

3 ( 1 

_|_ X I _ 1 _ 

2w k {eP u * — l)2cu g (e^ w « — 1) y (fc + q + p) 2 + m 2 

1 1 1 \ 

(k — q + p) 2 + m 2 (q — k + p) 2 + m 2 (k + q — p) 2 + m 2 J 

The second diagram in the r.h.s. of (6) is renormalized in the 2nd order of PT with the 
help of eqs. (6)-(9) i.e. 

\ = A + — _J_ f dAq 
0 2 (2tt) 4 7 (g 2 + m 2 ) 2 ’ 

2 o A 1 r d A q 

m ° = m ~ 2 (2^)4 J q 2 + m 2 ' 

Therefore np to the power A 2 we have 


m n — rn - 

U Q / 


^0 rn \ ' 1 


= —- tY, 
2 ^ 


\ ' _ 1 _ 

! I ,™2 _ o 1 2—/ (O^' 


(27t) 3 J q 2 + ml 2 ^ (27r) 3 J q 2 + m 2 

1 r d 3 q \ ( 1 f d 4 A; 




(27t) 3 J (q 2 + m 2 ) 2 ) \(2ir) 4 J k 2 + m 2 


3A 2 / 1 r d A q 
4 \(2n) 4 J (q 2 + rn 2 ) 2 


^ (27t) 3 J k 2 + m 2 



1 f d 3 q 

4 V ^ (27t) 3 J ( q 2 + m 2 ) 2 


rp y^ 1 / u ^ 

^ (27r) 3 J k 2 + m 2 


6 



e- 


A 2 2 [ d 3 qd 3 k 

6 ~^ q (27r) 6 J (g 2 + m 2 )(k 2 + m 2 )[(q + k — p) 2 + m 2 


(17) 


Ag 


Let us write the temperature dependent part of T®-function (4) as a sum of two terms 

7 (PI,P 2 ,T) = 72 (pl,p 2 ,T) + l3 (plp 2 ,T). 


Then combining (15)-(17) and regarding transformations from sums to integrals (10)-(14) 
we obtain 


72<yP 2 0 ,P 2 ,T) 


1 r d 3 q f A 2 1 [ d A k 

(27r) 3 J 2 u q {e^t — 1) { ^2 (27t) 4 J k 2 + m 2 * 


x 


(k — q — p) 2 + m 2 (k + q — p) 2 + m 2 k 2 + m 2 


where g 0 = 


(18) 


7 M,P 2 ,T) = 


d 3 qd 3 k 


2 !( 2 tt ) 6 J 2 UgieP”* - 1 ) 2 u k (e^ - 1 ) 


A 2 


(k + q + p) 2 + m 2 


+ 


+ 


(A; + g — p) 2 + m 2 (A; + p —g) 2 + m 2 (g + p — A;) 2 + m 2 q 


(19) 


where k 0 = iuOki % = 

One can recognize in the braces of (18) the renormalized 2^2 scattering amplitude at 
zero angle in the 2nd order of PT. So expected expression (5) is confirmed in 2nd order 
of PT. 

Now we assume by analogy that 73 should have the form 


( 2 -»2 rp\ _ 1 /' d 3 qd 3 kA^(p,q,k;p,q,k) 

7APo,P , ) 2 !(2tt ) 6 J 2u ^ q _ i)2ou k (e^ - 1) 


( 20 ) 


Really, the lowest order in which the 3^3 scattering amplitude appears is A 2 . There are 
10 tree diagrams (the direct channel diagram plus 9 exchanged diagrams) which contribute 
to the 3 —> 3 amplitude. Collecting them in 3 subsets as shown below 
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one can see that the first four terms in square brackets of (19) are just the four diagrams 
(21) with the allowance p — p', q — q', k — k! (zero angles scattering). More problematic 
are the diagrams (22). On mass shell of momenta q or k these diagrams diverge themselves 
and should be regularized. Several papers were devoted to this problem [13]. Our direct 
field-theory calculations give the simple recipe for such regularization. Namely one must 
set p — p' = e = (e 0 , e) and take the limit in the final result in the following order 


13(pI,P 2 ,T) = lim lim 


r=o ^ 0=0 2!(2vr) 3 


d 3 q d 3 k A^ (p, q, k; p + e, q, k )' 

2u q (ed u i — l)2uj k (e /3uJ k _ i) 


As to the last two diagrams (23), they do not contribute to y 3 because the considered 
r (2 kfunction by definition should contain 1P1 diagrams only. 


5 Conclusion 


Surprisingly the temperature dependent part of T^ 2 \pl,p 2 ,T) - eqs. (5) and (20) - can 
be expressed in an extremely compact and physically transparent form through the scat¬ 
tering amplitudes. These low order PT calculations give us the opportunity to generalize 
naturally enough the obtained result in the following way 

OO 

^(P 2 0 ,P 2 ,T) = J2~f n (p 2 0 ,p 2 ,T), 

71=2 


n\(2n] 


3 n 


n 


1n+l{pl,P 2 ,T ) = 

d 3 qi x 


(24) 


A (n+1 \p, q 1 ,...,q n ; p, q l ,...,q n ). 


2 - 1 ), 

This expression is linear in the scattering amplitudes and looks like a contribution of so 
called “first-type” diagrams in the S'-matrix formulation of statistical mechanics. We can 



not exclude without further analysis the appearence of some bilinear terms in amplitudes 
in the exact formula for 7(Po,p 2 , T) similar to the so-called “second-type” diagrams. But 
we shall not be surprised if the eq. (24) satisfies the exact representation. 

We thank profs. F. Paccanoni, B.V. Struminsky and E.S.Martynov for discussions. 
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